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Resistivity saturation is observed in many metallic systems with large resistivities, i.e., when the 
resistivity has reached a critical value, its further increase with temperature is substantially re- 
duced. This typically happens when the apparent mean free path is comparable to the interatomic 
separations - the Ioffe-Rcgcl condition. Recently, several exceptions to this rule have been found. 
Here, we review experimental results and early theories of resistivity saturation. We then describe 
more recent theoretical work, addressing cases both where the Ioffe-Regel condition is satisfied 
and where it is violated. In particular we show how the (semiclassical) Ioffe-Regel condition can 
be derived quantum-mechanically under certain assumptions about the system and why these 
assumptions are violated for high-T c cuprates and alkali-doped fullerides. 
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The electrical resistivity, p, of metals is usually calcu- 
lated in the Boltzmann theory, where the electrons are 
treated semiclassically. An electron is assumed to move 
with a wave vector k between scattering events, caused 
by phonons, other electrons, impurities or some other dis- 
order. The average distance an electron moves between 
two scattering events is the mean free path I. Assuming 
a three-dimensional system and a spherical Fermi surface 
with one sheet, p can be expressed in terms of / as (see 
Appendix A) 
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where Hf is the Fermi wave vector. Alternatively, if the 
resistivity is known experimentally, an apparent mean 
free path can be determined from Eq. (1). In the follow- 
ing, Eq. (1) will be used in the latter way to determine I, 
which then only depends on the experimental resistivity 
and the density of conduction electrons (giving fc^)- 
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FIG. 1. Resistivity of Cu, Nb 3 Sb (Fisk and Webb, 1976) and 
Nb (Abraham and Deviot, 1972). The figure also shows the 
the Ioffe-Regel (Ioffe and Regel, 1960) saturation resistivi- 
ties of NbsSb and Nb, setting the mean free path I in Eq. 
(1) equal to the distance between the Nb atoms. The cor- 
responding value for Cu, 260 ^f!cm, falls outside the figure. 
The figure illustrates that for NbsSb and Nb the resistivity 
saturates roughly as predicted by the Ioffe-Regel criterion, 
while p(T) ~ T for Cu at large T. 

For a good metal like Cu, I is of the order of hun- 
dreds or even thousands of A. In the semiclassical Boltz- 
mann theory and for T larger than a fraction of a typical 
phonon energy, the resistivity due to phonon scattering 
grows linearly with T. This is shown for Cu in Fig. 1. 
The apparent mean free path is then reduced correspond- 
ingly. For a good metal, however, I is much larger than 
the separation d of the atoms, even at the melting tem- 
perature. 

For certain metals, in particular many transition met- 
als and transition metal compounds, the resistivity be- 
haves in a completely different way. This was emphasized 
by Fisk and Webb (1976), who studied the A15 com- 
pounds Nb 3 Sb (see Fig. 1). Similar results had earlier 
been obtained by Woodard and Cody (1964) for NbsSn 
and by several other groups for other compounds. For 
small T, p(T) grows much faster than for Cu and I be- 
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comes comparable to d already for T of the order of sev- 
eral hundred K. At this point, the slope of the resistivity 
curve is substantially reduced, and I stays comparable 
to d for experimentally accessible values of T. This is 
referred to as resistivity saturation. It describes the sit- 
uation where p(T) grows much slower than p(T) <~ T, 
predicted by the Boltzmann equation, but it does not 
necessarily mean that p{T) becomes a constant. 

Ioffe and Regel (1960) pointed out that the semiclassi- 
cal theory makes no sense if I < d, and we refer to I > d 
as the Ioffe- Regel condition. Inserting I — d in Eq. (1) 
gives the Ioffe- Regel resistivity, shown in Fig. 1. In a 
semiclassical picture, it is natural that I cannot be much 
smaller than d, since one may expect that an electron 
at most is scattered at every atom. Saturation is then 
expected when I ~ d. A semiclassical theory, however, 
breaks down when I ~ d, since the uncertainty in the 
k-vector of an electron is comparable to the size of the 
Brillouin zone. A semiclassical theory cannot therefore 
explain saturation. Nevertheless, as discussed in Sec. II, 
in the 1970's and early 1980's a large number of metals 
were found which showed saturation when I ~ d, and the 
behavior seemed to be universal. During this time, much 
theoretical work was performed to explain saturation. No 
consensus was reached, however, and the interest turned 
to other problems. 
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FIG. 2. Resistivity of Bi 2 Sr 2 Cai- :c Y x Cu208+ ! , (T c = 30 
K) (Wang et al, 1996ab), Lai.93Sro.o7Cu04 (Takagi et al, 
1992), Ndi.84Ce .i6Cu4-» (T c = 22.5 K) (Hikada and Suzuki, 
1989), YBa 2 Cu306 +;c (T c = 60 K) (Orenstein et al, 1990), 
Bi 2 Sr 2 Cu 6+! , {T c = 6.5 K) (Martin et al, 1990) and Nb 3 Sb 
(Fisk and Webb, 1976). The arrow shows the Ioffe- Regel re- 
sistivity of Lai.93Sro.o7Cu04- The figure illustrates that there 
is no sign of saturation at the loffe-Regel resistivity, but in 
some cases perhaps at much larger resistivities. Observe the 
magnitude compared with NbaSb. 

The situation changed drastically in the late 1980's, 
when it was found that the high-T c cuprates behave very 
differently (Gurvitch and Fiory, 1987). Fig. 2 shows 



that the resistivity for these compounds is typically much 
larger than for metals as NbsSb, satisfying the loffe-Regel 
condition. As described in Appendix A, the loffe-Regel 
resistivity for L^-zSr^CuO^ is very large, about 0.7 
mficm, due to the small carrier concentration. Never- 
theless, the resistivities of the high-T c cuprates substan- 
tially exceed the Ioffc-Rcgcl resistivity or are compara- 
ble to this resistivity without signs of saturation. The 
Ioffc-Rcgcl condition is therefore strongly violated. Both 
La2- :E Sr :I ;Cu04 for small x and Bi2Sr2Cai_ :E Y ;r Cu208+j, 
show, however, signs of saturation, although at much 
larger values than the loffe-Regel resistivity. Fig. 3 shows 
other examples of the violation of the Ioffc-Rcgel condi- 
tion, namely Rb 3 Cgo (Hebard et al, 1993), La 4 Ru 6 0ig 
(Khalifah et al, 2001) and Sr 2 Ru0 4 (Tyler et al, 1998). 
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FIG. 3. Resistivity of Rb 3 C 6 o (Hebard et al, 1993) 
La 4 Ru 6 Oi 9 (Khalifah et al, 2001), Sr 2 Ru0 4 (Tyler et al, 
1998), Nb 3 Sb (Fisk and Webb, 1976) and the loffe-Regel re- 
sistivity for RbsCeo- There is no sign of saturation at the 
Ioffe- Regel resistivity, but La 4 Ru 6 0ig may saturate at a much 
larger resistivity. 



Although there are substantial uncertainties in the ab- 
solute values of the experimental resistivity, it is, never- 
theless, clear that the loffe-Regel condition can be vio- 
lated. This shows that the semiclassical argument for sat- 
uration is not just questionable, but actually gives wrong 
predictions for these systems. This emphasizes the need 
for a theory of saturation going beyond the semiclassical 
treatment of the electrons. Such a theory should also 
explain the violation of the loffe-Regel condition for the 
systems mentioned above. The drastic change in the ex- 
perimental situation over the last 15 years has led to a 
renewed interest in resistivity saturation. 

In Sec. II, we review the experimental results, and in 
Sec. Ill early theoretical work. We then briefly describe 
methods for calculating the resistivity in Sec. IV. Much 
of the theoretical work has started from a Boltzmann- 
like approach, assuming that the scattering mechanism 
is a relatively weak perturbation. In Sec. V, we de- 



2 



scribe an approach starting from the opposite limit of 
very strong scattering, assuming that the Drude peak 
has been completely removed. This approach is based 
on the f-sum rule. It leads to an approximate upper 
limit /9 sa t to the resistivity, which usually has a weak 
T dependence. Resistivity saturation then happens if 
the initial slope of p(T) is so large that p sa ,t is reached 
at small values of T and if p sa t has a weak T depen- 
dence. This approach is applied to three models with 
different saturation behavior. In Sec. VI, we describe 
the treatment of a model of weakly correlated electrons 
in a broad band, appropriate for many transition metal 
compounds. This model shows saturation in agreement 
with the Ioffe-Regel condition. We show how this condi- 
tion can be derived quantum-mechanically by assuming 
i) noninteracting electrons and ii) T -C W, where W is 
the band width. In Sec. VII we treat high-T c cuprates, 
for which assumption i) about noninteracting electrons is 
violated. Using the t — J model, we find that saturation 
can occur, but at much larger resistivities than predicted 
by the Ioffe-Regel condition. Other metals violating the 
Ioffe-Regel condition are discussed in Sec. VIII. In Sec. 
IX, we describe a model of alkali-doped Cgo , which shows 
no saturation. This result depends on the intramolecular 
character of the phonons and on condition ii) (T <C W) 
being violated. The relation to Anderson localization and 
Mott's minimum conductivity is discussed in Sec. X. 
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FIG. 4. Resistivity of Tii-^Al^ alloys. The figure suggests 
that the saturation resistivity is independent of the disorder 
(after Mooij, 1973). 



II. EXPERIMENTAL RESULTS 

Most cases of resistivity saturation have been found for 
transition metal compounds. Mooij (1973) made one of 
the first observations for Tii^Al^ (see Fig. 4) and sev- 
eral other alloys. In addition, he found that the positive 
T dependence of the resistivity becomes weaker or even 
negative for strongly disordered systems. Similar results 
were found by Arko et al. (1973) and by Tsuei (1986) 
for other compounds. These results suggest that compo- 



sitional disorder and disorder due to thermally excited 
phonons can have a similar effect. Other early examples 
were found by Fisk and Lawson (1973). 

As discussed in the introduction, resistivity saturation 
was found for A15 compounds (Woodard and Cody, 1964; 
Marchenko, 1973; Fisk and Webb, 1976; Wiesmann et 
at, 1977; Gurvitch et at, 1978). Another class of sys- 
tems showing resistivity saturation is the Chevrel phases 
(Martin et al., 1978; Sunandana, 1979). Early measure- 
ments gave very large values for the resistivity, probably 
because of sample problems, but later measurements on 
single crystals gave resistivities comparable to those of 
A15 compounds (Kitazawa et at, 1981). Fig. 5 shows the 
resistivity of the 3d and 5d transition metals. Several of 
these metals show clear signs of saturation, in particular 
a-Mn. Other cases are less clear, and, e.g., the resistivity 
of W even bends slightly upwards, although it is already 
large. As discussed at the end of this section, W is nev- 
ertheless an example of saturation, if anharmonic effects 
are considered. The above metals can be considered as 
cases where the Ioffe-Regel condition is satisfied. 




T(K) 



b) 


Hi/" 










Ta 




.X' Re ... .. 








W 


La// 














" : V V :> Pt 






Os 











500 1000 1500 2000 2500 3000 3500 
T(K) 

FIG. 5. Resistivity of 3d (a) and 5d (b) transition metals 
(Bass, 1982). 

We next consider the in-plane resistivity of high-T c 
cuprates. The resistivity depends strongly on the dop- 
ing of the Cu02 plane. Information about the doping 
is given by the superconductivity transition temperature 
T c , which is quoted below, when available. As the dop- 
ing is reduced below the optimum value, T c drops and p 
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increases. In the case of I^-rrSr^CuC^ we assume that 
the doping is given by x. Results for several cuprates 
were given in Fig. 2, and we give some further exam- 
ples in Table I. These are examples of systems where the 
resistivity is comparable to the Ioffe-Regel resistivity al- 
ready at moderate values of T and without showing signs 
of saturation. 

TABLE I. Resistivity p(T) (in mQcm) of high-T c cuprates. 
The measurment temperature T and the superconductivity 
transition temperature T c are given in K. 



Compound 



T c T p(T) Reference 



HgBa 2 CaoCui04+^ 94 

HgBa 2 CaiCu 2 06 + ^ 122 

HgBa 2 Ca 2 Cu 3 8+:c 125 

HgBaaCasOuOio+s 130 

Tl 2 Ba 2 Cu0 6+!/ 80 

Tl 2 Ba 2 Cu0 6+!/ 80 

TlSr 2 CaCu 2 7 - y 65 

Bi 2 Sr 2 CaCu208 +y 76 



300 0.5 Daignere et al, 2001 

300 0.3 Yan et al, 1998 

500 0.6 Carrington et al, 1994 

400 0.5 Lohle et al, 1996 

300 1.3 Kubo et al, 1991 

270 0.6 Duan et al, 1991 

300 0.5 Kubo et al, 1991 

300 1.2 Chen et al, 1998 



Eh^S^Cai-zY^C^Os+y is of particular interest, since 
Wang et al. (1996ab) found saturation for the most un- 
derdoped samples, e.g., x = 0.11 (T c = 30 K) as shown 
in Fig. 2. Similar results have also been obtained by 
Forro et al (2002). Other measurements only showed 
weak signs of saturation (Mandrus et al., 1992; Ruan et 
al, 2001). In the measurements of Ruan et al. (2001), 
the reason may have been that their x = 0.11 sample had 
a higher T c = 56 K than that of Wang et al. (1996ab) 
(T c = 30 K), suggesting that the former sample was less 
underdoped. We conclude that the Ioffe-Regel condition 
is violated for high-T c superconductors but that some 
systems with very large resistivities may show signs of 
saturation at much larger values than the Ioffe-Regel re- 
sistivity. 
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FIG. 6. Resistivity of CaRuOs (Klein et al, 1999ab), Cr0 2 
(Rodbell et al, 1966), V0 2 (Allen et al, 1993) and SrRuOs 
(Allen et al, 1996). 

Apart from the high-T c cuprates, a number of tran- 



sition metal compounds have been found which violate 
the Ioffe-Regel condition. Fig. 3 and Fig. 6 shows some 
examples. Some of these metals show signs of satura- 
tion (La4RugOig, Cr02 and perhaps CaRuOs), but at 
higher resistivities than the Ioffe-Regel resistivity, while 
other systems show no sign of saturation. Several of the- 
ses metals are believed to be non- Fermi liquids. The 
data in Figs. 3 and 6 were obtained from single crys- 
tals, except for RbjCeo and CaRuC>3. This could be a 
reason why the resistivity in Fig. 6 is much larger for 
CaRuOs than for the related SrRuC>3. Large resistivities 
are also found in manganites (Salamon and Jaime, 2001). 
These systems raise interesting problems in terms of the 
colossal magnetoresistance, possibly important polaron 
effects and many phase transitions, but they are outside 
the scope of this paper. 

The resistivity of A 3 Ceo (A= K, Rb) shows a quite sub- 
stantial spread between different measurements (Hebard 
et al, 1993; Hou et al, 1993; Palstra et al, 1994; Degiorgi 
et al, 1994). These have been performed for thin films 
and doped single crystals using direct and optical meth- 
ods. In particular the thin films may have substantial 
defects (e.g., grain boundaries) which would increase the 
resistivity. The optical measurements should, however, 
be much less sensitive to this. Indeed, optical measure- 
ments for doped single crystals (Degiorgi et al, 1994) 
gave resistivities that are about a factor of two smaller 
than those obtained from direct measurements for thin 
films (Palstra et al, 1994), but comparable to (some- 
what larger than) direct measurements for doped single 
crystals (Hou at al., 1993). This suggests that the thin 
film results in Fig. 3 could overestimate the resistivity of 
Rb 3 C6o by a factor of two or somewhat more. This would 
not change the qualitative conclusion that the Ioffe-Regel 
condition is violated for A^Cqo. 

Experiments for Rb 3 C6o show no signs of saturation 
up to about 500 K (Hebard et al, 1993). It is difficult to 
reach higher values of T, due to the possibility of ther- 
mally driven rearrangements of the alkali atoms. Hou et 
al (1995), using a pulsed heating technique to reach 800 
K, found a small change in the slope at about T = 500 K. 
This was interpreted as a sign of saturation. Using the 
parallel resistor formula (Wiesmann et al, 1977), they 
deduced the saturation resistivity 6 ± 3 mfkm, corre- 
sponding to the mean free path I = 1 ± 0.5 A. Thus there 
is no saturation at I ~ d, where d = 10 A is the separation 
of two Ceo molecules (see Appendix A) , but it is hard to 
judge if there is saturation at some larger resistivity. 

The Boltzmann equation predicts p(T) ~ T for large T 
only if various parameters, such as the electron-phonon 
interaction, are independent of T. This assumption is 
sometimes strongly violated. For instance, anharmonic 
effects reduce the phonon frequencies of W by about 30 
% at the melting point (Grimvall, 2001; Grimvall et al, 
1987; Guillermet and Grimvall, 1991). According to the 
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Boltzmann equation, this should increase p(T) by about 
a factor of two, due to an increase in the electron-phonon 
coupling, giving a growth of p(T) which is much faster 
than linear. Experimentally, however, the growth is only 
slightly faster than linear (Fig. 5b). This suggests that 
saturation in W is masked by the reduction of the phonon 
frequencies (Grimvall, 2001). 

The resistivity is usually measured at constant pres- 
sure. As T is increased, the crystal expands and vari- 
ous parameters arc changed. Sundqvist and Andersson 
(1990) converted the resistivity of YBa2Cu307_ y to con- 
stant volume data, which is the more relevant quantity 
for theoretical discussions. They found that this lowered 
the T — 500 K resistivity by about 32 %. Although 
this conversion involves many uncertainties (Sundqvist 
and Andersson, 1990), it was concluded that the pro- 
nounced linearity in the constant pressure data is lost. 
Actually, the corrected data show signs of saturation. A 
similar correction of the data for RbaCgo (Vareka and 
Zettl, 1994) changed the approximately quadratic T de- 
pendence found for constant pressure (see Fig. 3) to an 
approximately linear dependence for constant volume. 

III. EARLY THEORETICAL WORK 

Resistivity saturation was studied theoretically very 
intensively from the middle of the 1970's to the early 
1980's. In this section we review some of this work 
as well as some later work. We focus on work related 
to the electron-phonon scattering, while work based on 
electron-electron scattering (Jarrell and Pruschke, 1994; 
Lange and Kotliar, 1999; Parcollet and Georges, 1999; 
Merino and McKenzie, 2000; Gunnarsson and Han, 2000) 
fall outside the scope of this review. 

Much of the early work has been reviewed by Allen 
(1980ab). Initially, explanations were proposed in terms 
of strong variations of the electronic structure on the en- 
ergy scale of ksT or unusual anharmonic effects. These 
effects appear, however, to be too small to explain exper- 
iment, and as ever more examples of resistivity satura- 
tion were discovered, such theories had to be abandoned 
(Allen, 1980ab). 

Wiesmann et al. (1977) found empirically that the 
resistivity of systems such as NbaSb can be rather accu- 
rately described by a parallel resistor formula 



P(T) Psat Pideal(T)' 

where p sat is a saturation resistivity and Pidea,i(T) is the 
resistivity of the Boltzmann equation, i.e., linear in T for 
large T. 

Mooij (1973) observed that strong disorder can lead 
to a negative slope of p(T) (Fig. 4). Jonson and Girvin 
(1979), Girvin and Jonson (1980) and Imry (1980) ar- 
gued that this is due to an incipient Anderson transition. 



The disorder is assumed to be so large that the system 
is close to an Anderson transition for T = 0. As T is 
increased, inelastic scattering (phonon assisted hopping) 
leads to a loss of phase information (Lee and Ramakrish- 
nan, 1985) and the system moves away from the Ander- 
son transition. This was argued to lead to a reduction of 
the resistivity, as observed by Mooij (1973) for strongly 
disordered systems. 

Weger and Mott (1985) and Weger (1985) proposed a 
theory in terms dchybridization of the s- and d-electrons 
at large values of T. Laughlin (1982) presented a theory 
where exchange interaction reduced the density of states 
and increased the Fermi velocity in such a way that re- 
sistivity saturation was obtained. 

Cote and Meisel (1978) and Morton et al. (1978) pro- 
posed that an electron is only scattered by a phonon if the 
mean free path of the electron is longer than the phonon 
wave length. As T is increased and I becomes shorter, an 
increasing number of phonons become inefficient scatters. 
With these assumptions, resistivity saturation can be de- 
rived. Later calculations treating the electrons quantum 
mechanically, as described before Eq. (5), do not seem 
to support the basic assumption in this theory (Calandra 
and Gunnarsson, 2002). 

Several theories have been based on models with one 
electronic state per unit cell. Christoph and Shillcr 
(1984) studied the Frohlich Hamiltonian. Using a trun- 
cated equation of motion approach, they included some 
higher order terms in the electron-phonon interaction. In 
this way they derived the parallel resistor formula (Eq. 
(2)). As pointed out by them, their p sat ~ 0.84 mficm is 
substantially larger than what is observed for transition 
metal compounds, but it is comparable to what might 
be expected for a one band model, as in Eq. (11) be- 
low (Nd = 1). Bclitz and Schirmacher (1983) studied a 
generalized Frolich model which also included impurity 
scattering. They included higher order terms by using a 
mode coupling technique in a memory function approach. 
They found a rather good agreement with experiment, 
but a substantially smaller p sat than might have been 
expected from Eq. (11) for = 1. 

Ron et al. (1986) treated a Kronig-Penney chain, 
where the positions of the (5-functions vibrate. Using the 
Landaucr formula they found saturation in the limit of 
large amplitude vibrations. Their p sat depends, however, 
on the strength of the potentials, and it therefore appears 
that the Ioffe-Regel condition is not explained. 

The Boltzmann equation only includes intraband tran- 
sitions. Chakraborty and Allen (1979) and Allen and 
Chakraborty (1981) emphasized that resistivity satura- 
tion is typically found for systems where interband tran- 
sitions are important. They included these transitions 
by generalizing the Boltzmann equation, treating terms 
of the order (d/l)°, while the Boltzmann equation only 
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contains terms of the of the order (d/l) -1 . They demon- 
strated that this leads to a parallel resistor type of for- 
mula (Eq. (2)) and estimated the saturation resistivity to 
be of the right order of magnitude. Due to the complex- 
ity of the equations, it was not possible to make explicit 
calculations and the sign of the p sat in Eq. (2) could not 
be determined. They also pointed out that higher order 
terms in (d/l) may be important. 
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FIG. 7. Resistivity of the C60 model for different electron- 
phonon couplings A according to QMC calculations (Gunnars- 
son and Han, 2000; Calandra and Gunnarsson, 2002). The 
phonon frequency is u) v h = 0.1 eV. The figure illustrates the 
lack of saturation. For A = 0.80 the onset of superconductiv- 
ity can be seen as a sharp downturn in p(T) as T is lowered, 
due to superconducting fluctuations. For A = 1.06 and 1.32, 
the resistivity has a negative slope for small T, indicating an 
insulating system. 



Millis et al. (1999) studied a Holstein like model, where 
the phonons couple to the electron level energies (in the 
following called LE coupling). They found that there 
is a change of slope in the p{T) curve for the case of 
a strong electron-phonon coupling, and associated this 
with resistivity saturation. They observed, however, that 
the resistivity does not saturate and remarked that resis- 
tivity saturation is a misnomer. Gunnarsson and Han 
(2000) performed a quantum Monte-Carlo (QMC) calcu- 
lation for a somewhat related model of A 3 Ceo (A= K, 
Rb). This model included the partly occupied, three- 
fold degenerate t\ u level as well as the LE coupling to a 
five-fold degenerate H g phonon. We refer to this as the 
Ceo model. The results, shown in Fig. 7, are somewhat 
similar to those of Millis et al. (1999), but they were con- 
sidered as an example of lack of saturation. It was later 
shown (Calandra and Gunnarsson, 2001) that a more re- 
alistic description of saturation can be obtained by using 
a coupling to the hopping integrals (HI coupling) (see 
Sec. VI). 



IV. CALCULATION OF RESISTIVITY 

The resistivity is usually calculated in the Boltzmann 
theory. The electrons are treated semiclassically, and are 
assumed to be accelerated by the electric field between 
scattering processes, caused by impurities, phonons or 
other electrons. This can be considered as the lowest 
order theory in an expansion in l/(kpl) (Kohn and Lut- 
tinger, 1957). Both because of this and because of its 
prediction p(T) <~ T for large T, it is clear that the Boltz- 
mann equation cannot be used for describing saturation. 
Instead, it is necessary to treat the electrons quantum- 
mechanically. This can be done in the Kubo formalism. 
For an isotropic system, this requires the calculation of 
a current-current correlation function (Mahan, 1990) 



1 [*> 



dre™ (T T j(<r) ■](<))>, 



(3) 



where N is the number of atoms, ft is the volume per 
atom, (3 — l/(fcsT), kg is the Boltzmann constant, T T 
is a time-ordering operator and j is the current operator. 
ir(iu) is analytically continued to real frequencies, giving 
7r rot (w). The optical conductivity is then 



a(u>) 



Im7r rot (w) 



(4) 



and the resistivity is l/cr(0). 

ir(iu>) can be calculated for imaginary times by using 
a determinantal quantum Monte-Carlo (QMC) method 
(Blankenbcclcr et al, 1981). It is then analytically con- 
tinued to real frequencies by using a maximum entropy 
method (Jarrell and Gubernatis, 1996). This approach is 
very useful for establishing the saturation properties of a 
given model. 

To interpret the results, a simpler method is used, 
where the phonons are treated semiclassically (Calandra 
and Gunnarsson, 2001; 2002). A large supercell with N 
atoms and periodic boundary conditions is considered. 
Each phonon coordinate is given a random displacement 
according to a Gaussian distribution with the width de- 
termined by T. This set of displaced coordinates (snap 
shot) defines a new Hamiltonian. While the original 
Hamiltonian represents a complicated many-body prob- 
lem due to the electron-phonon interaction, the dynamics 
of the phonons have been eliminated in the new Hamilto- 
nian. Since we assume no electron-electron interaction, 
this is then a one-particle Hamiltonian for a disordered 
system. The corresponding eigenstates \l) and eigenval- 
ues e; are found. The Kubo-Greenwood formula (Green- 
wood, 1958) is then used 



2tt 



J2\(i\h\i)\ 2 (fi-MS(tuv 



u' 



(5) 



G 



where j x is the current operator in the x direction, and 
/; is the Fermi function for the energy It is essential 
that in this approach only the phonons are treated semi- 
classically and that the electrons are treated quantum- 
mechanically. This is in contrast to the semiclassical 
treatment of the electrons in the Boltzmann theory, 
which is unable to describe saturation. 

In the following sections we discuss a new theory for 
describing different classes of materials with different sat- 
uration properties (Gunnarsson and Han, 2000; Calandra 
and Gunnarsson, 2001; 2002; 2003). 
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FIG. 8. o(u>) as a function of ui for the NbsSb model in the 
semiclassical treatment of the phonons. The frequency has 
been scaled by the T = band width W. The figure illus- 
trates how the Drude peak is lost as T is increased. 



V. f-SUM RULE 

For small T, the optical conductivity cr(u>) typically 
has a Drude peak, due to intraband transitions between 
quasiparticles and centered at u> — 0, and a remaining 
incoherent part, as illustrated in Fig. 8 for NbsSb. As T 
is increased, the height of the Drude peak is reduced, and 
for large T it disappears. This is also illustrated in Fig. 
9, which shows the experimental optical conductivity of 
Cr02 at T = 300 K. The arrow marks the conductivity 
corresponding to the saturation resistivity (see Fig. 6). 
This conductivity is similar to the magnitude of the inco- 
herent contribution (at 300 K). Again, this suggests that 
as T is increased, the Drude peak disappears without the 
incoherent part changing very much. In the following, we 
therefore assume that T is so large that the Drude peak 
has disappeared (Calandra and Gunnarsson, 2002). 

The Boltzmann theory starts from the periodic system 
and treats the scattering mechanisms as small perturba- 
tions. The theory focuses on intraband transitions and 
the Drude peak. Here we start from the the opposite 
limit, where the scattering is so strong that the Drude 
peak disappears. 

To estimate the incoherent contribution to <j(lo), we 
use the f-sum rule. We focus on tight-binding models 



with one type of electrons, e.g., d electrons in a transi- 
tion metal compound. Then there are no on-site matrix 
elements of the current operator. The off-site matrix el- 
ements can be related to hopping matrix elements, using 
charge and current conservation. For an isotropic three- 
dimensional system, the tight-binding version of the f- 
sum rule then takes the form (Maldague, 1977; Calandra 
and Gunnarsson, 2002) 

where Ek is the hopping (kinetic) energy, N is the num- 
ber of atoms and Q is the volume per atom. For a system 
with just nearest neighbor hopping and hxed atoms, d is 
the separation of the sites. If the atoms vibrate or if 
there is hopping between more distant neighbors, d is a 
weighted average of neighboring distances (Calandra and 
Gunnarsson, 2002). 




co(eV) 

FIG. 9. <j(lu) for Cr02 according to experiment (full line) at 
T = 300 K (Singley et al, 1999) and the large T model in Fig. 
10 (dotted line), adjusting the area to the sum rule (6). The 
experimental intensity above oj ~ 1.5 eV is due to interband 
transitions. The arrow shows the conductivity corresponding 
to the saturation resistivity in Fig. 6. It suggests that for 
large T, the Drude peak in the experimental a(u) disappears. 

For models with only one type of electrons, e.g, d elec- 
trons, <t(lu) at large T looks schematically as in Fig. 10. 
This can be compared with the calculation for NbsSb in 
Fig. 8 or the expected large T contribution from the ti g 
band of Cr02 in Fig. 9 (Drude peak gone), excluding 
the interband transitions for fou) > 1.5 eV. For ftu larger 
than the band width W, <t(oj) 0. If <j(oj) furthermore 
lacks large structures around uj = 0, we can approximate 

a[uj) ~ \ 0, otherwise. (?) 

Then J a(uj)hdui = a(0)W/j, and the f-sum rule (6) re- 
quires that <r(0) is given by the right hand side of Eq. (6) 
multiplied by jh/W (Calandra and Gunnarsson, 2002; 
2003), i.e., 
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a(0) = 



w 



<j((jj)du 



2„2 



7T7 d z e 

ewmrn 



\Ek\, 



(8) 



where 7 ~ 2 depends on the shape of a(ui). This provides 
an approximate lower limit to the conductivity (upper 
limit to the resistivity) , which may be reached if T is so 
large that the Drude peak is negligible. The conductivity 
could be appreciably smaller only if there is a substantial 
amount of weight for lo > W or if <j(u>) is anomalously 
small for ui = 0. The former should not be important, 
since Hio > W requires multiple electron-hole pair excita- 
tions. Since the current operator is a one-particle oper- 
ator, these excitations should have a small weight. The 
latter could happen in the case of an incipient Anderson 
transition, but it should normally not occur at large T, 
where inelastic scattering destroys the phase information 
required for an Anderson localization (Lee and Ramakr- 
ishnan, 1985). The (lower) limit (8) to the conductivity 
is T dependent, due to the T dependence of Ek and W. 
This T dependence is often, but not always, rather weak. 
Below we use Eq. (8) to discuss different classes of met- 
als. 



o(0) 



a(0)/y 




FIG. 10. Schematic picture of a(u>). The average of a(ui) over 
the band width is given by a(0)/j, where 7 ~ 2. 



VI. WEAKLY CORRELATED BROAD BAND SYSTEMS 

We first focus on metals with a large band width W. 
We therefore assume that i) W is so large that the elec- 
trons can be treated as noninteracting and ii) T <C W . 
In particular, we focus on NbaSb, showing a pronounced 
saturation, and Nb, showing a much weaker saturation 
(see Fig. 1). We use a model with an Ad-fold degenerate 
orbital on each site. The orbitals on different atoms are 
connected by hopping matrix elements, t im - m ', where i 

and j label the atoms and m,m' = 1, ..,Nd the orbital. 
Then the electronic Hamiltonian is 



HCl = E ^, 

i^j,mm a 



jm vma jm <r' 



(9) 



Three Einstein phonons are introduced for each atom, 
describing its displacement in the three coordinate di- 
rections. The phonons couple to the hopping integrals, 



since these depend on the atomic positions. We refer to 
this as HI coupling and the model is referred to as the 
transition metal (TM) model. The model for Nb has 4c? 
orbitals (Nd = 5) on a bec lattice. For NbaSb we use the 
NbJ model, neglecting the Sb atoms (Pickett et at, 1979), 
and putting Ad orbitals on the Nb sites. The models are 
therefore identical except for the lattice structure. For 
strongly ionic systems, the coupling to the level energies 
(LE) may also be important. 

Fig. 11 shows the QMC (circles) results for the resis- 
tivity of Nbg model for a supercell with N = 36 atoms. 
These results extrapolate to a large resistivity at T = 0. 
On the other hand, the resistivity of the model must be 
zero at T = 0. It then follows that there must be a 
large change of the slope of p(T) as T is increased, i.e., 
resistivity saturation. This model can therefore be used 
for analyzing saturation. Fig. 11 also shows that the 
semiclassical treatment of the phonons is quite accurate 
for the TM model, and that it therefore can be used for 
interpreting the results. 

To use f-sum rule of Sec. V, we calculate the hopping 
energy Ek- Since we have assumed that the electrons 
are noninteracting and that T <C W, we obtain 



E K = 2N d / eN{e)de = -2aN d WN, (10) 

J-W/2 

where Ep is the Fermi energy, N(e) is the density of 
states per atom, orbital and spin, A^ is the orbital degen- 
eracy, and a depends on the shape of N(e) and the filling. 
Typically a <~ 0.1 for a metal close to half-filling. The 
dependence on filling is moderate around half-filling (Ca- 
landra and Gunnarsson, 2002), but a = for an empty 
or full band. Inserting Eq. (10) in Eq. (8) gives an ap- 
proximate lower limit to the conductivity (Calandra and 
Gunnarsson, 2002; 2003) 



o- S at(0) 



7ra7 d 3 N d e 2 

ITU hd 



(ii) 



Here 7ra7/3 <~ 0.2 depends on the details of the electronic 
structure, and d 3 /0 ~ 1 depends on the lattice structure. 
The result is independent of W, since the W in Eq. (8) 
cancels that in Eq. (10). The essential material param- 
eters, d and N d appear in N d e 2 / (Hd), which has the di- 
mension conductivity. The corresponding apparent mean 
free path is (Calandra and Gunnarsson, 2001) 



cNjd, 



(12) 



where c ~ 0.6 — 0.8, which is an approximate lower limit 
to /. It then follows that I > d. This provides a quantum- 
mechanical derivation of the Ioffe-Regel condition, as- 
suming that the electrons can be treated as noninteract- 
ing and that T -C W. For a transition metal compound, 
with N d = 5 and d <~ 3 A, this leads to a resistivity of 
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the order of 0.1-0.2 mficm, in agreement with the exper- 
iments. 
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FIG. 11. Resistivity for the Nb^ model of NbaSb, comparing 
the semiclassical [ broken (N — 36 atoms) and full (N = 648) 
curves] and QMC (circles, N = 36) treatment of the phonons. 
The figure also shows the small (Eq. (13)) and large (Eq. 
(11)) T results. It illustrates that the resistivity of the TM 
model saturates at large T (after Calandra and Gunnarsson, 
2001). 

We next consider the small T behavior of p(T). In this 
limit the Boltzmann theory can be used, since I 3> d. For 
T larger than some fraction of u p h (Grimvall, 1981) 

,XD = 8 ^. (13) 

where A is the dimensionless electron-phonon coupling 
constant and fl p i is the plasma frequency 

(^) 2 = ^E J B /k[^f6(e„W-E F ). (14) 

Here £„(k) is the energy of a state with the band index 
v and the wave vector k. O p ; depends on the average 
Fermi velocity. 

The straight line corresponding to Eq. (13) is shown 
in Fig. 1 1 . It shows how the conductivity is reduced and 
the resistivity is increased due to the reduction of the 
Drude peak height as T is increased. The horizontal line 
in Fig. 11 corresponds to the large T limit (Eq. (11)) 
of the conductivity. When this horizontal line is crossed 
by the steep line (Eq. (13)), the Drude peak has been so 
strongly reduced that it disappears under the incoherent 
part, and is spread out over a large energy range. As 
discussed above, due to the f-sum rule, the conductivity 
cannot normally be much further reduced at this point. 

Saturation requires that the two lines (Eqs. (11, 13)) 
cross in the available temperature range. For most met- 
als, the line in Eq. (13) has such a small slope that the 
crossing would happen far above the melting point. 

The model of Nb shows a much weaker saturation 
than the Nb^ model of NbaSb (Calandra and Gunnars- 
son, 2001), as is also found experimentally (see Fig. 1). 
This is due to a much larger O p ; for Nb (9.5 eV in the 



model) than for NbsSb (3.4 eV). The slope of the line 
in Eq. (13) is therefore a factor of five larger for Nb 3 Sb 
than for Nb. This leads to the much more pronounced 
saturation for NbjSb. The difference in f2 p ; is due to 
the large unit cell of Nbg, which leads to many bands 
and many forbidden crossings. This results in quite flat 
bands and small electron velocities and a small O p ; (Eq. 
(14)). In this context it is interesting to note that a-Mn 
has the most pronounced saturation among the transi- 
tion metals (see Fig. 5a). The unit cell of a-Mn has 
58 atoms, suggesting a very small tt p i and a very steep 
initial slope of p(T). Chakraborty and Allen (1979) and 
Allen and Chakraborty (1981) observed that saturation 
usually happens for systems with important interband 
transitions. Such systems typically have large unit cells 
and small fi p j . Their observation can then be understood 
in terms of the arguments above. 

For a one-band model with one s-state per unit cell, 
/3 sa t is large (Nd = 1 in Eq. (11) instead of — 5 for the 
TM model) and X/fl p i is typically small. The crossing of 
the two lines in Eqs. (11, 13) then happens for such large 
vibration amplitudes that a calculation is not very mean- 
ingful. The main reason is that for an s-band model, the 
matrix elements do not change sign as the atoms vibrate, 
and very large vibration amplitudes are needed to make 
the scattering so strong that the Drude peak essentially 
disappears. It is therefore interesting to study a one- 
band model where the state has many strongly directed 
lobes. Even moderate atomic vibrations then change the 
signs of the matrix elements. As a model we use one 
of the three t\ u orbitals of a Cgo molecule (hypothetical 
case of a strong crystal-field splitting). This orbital has 
I = 5 and correspondingly many nodes. The resulting 
resistivity saturation is moderate, with the slope at large 
T being reduced by about a factor of three to four in the 
semiclassical treatment of the phonons. This illustrates 
that interband transitions are not indispensable for resis- 
tivity saturation. 

The optical conductivity in Figs. 8 and 9 can be viewed 
as a Drude peak riding on top of a broad structure of in- 
coherent contributions. According to the f-sum rule and 
the arguments above, the height and shape of the latter 
part does not change much for values of T discussed here. 
We then approximately have 



<t(uj = 0, T) = er Drudo (w = 0,T) + cr sat , (15) 



where ODrudc(^ — 0, T) is the inverse of the resistivity in 
Eq. (13) and cr sa t is the height of the incoherent contri- 
bution (Eq. (11)). This is the parallel resistor formula 
(Wiesmann et al, 1977). 
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VII. STRONGLY CORRELATED SYSTEMS. HIGH-T C 
CUPRATES 

In the derivation of the Ioffe-Regel condition, we as- 
sumed i) noninteracting electrons. We now focus on the 
high-T c cuprates, where this assumption is invalid. The 
transport properties in the CuC>2 plane are mainly deter- 
mined by the antibonding band of Cu x 2 — y 2 and O 2p 
character. It is then natural to consider one (x 2 — y 2 ) or- 
bital per site and the orbital degeneracy = 1. These 
orbitals are put on a two-dimensional square lattice, de- 
scribing the CuC>2 plane. The Coulomb interaction is 
described by a Hubbard U acting between two electrons 
on the same site. U is assumed to be so large that states 
with double occupancy on a site can be projected out, 
leading to the t — J model (Zhang and Rice, 1988). 
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FIG. 12. In-plane optical conductivity of L^-xSr^CuO^ as 
a function of u) for x = 0.06 and x = 0.10 (Takenaka et al, 
2002). The figure illustrates how the peak at u> — gradually 
goes away as T is increased and that the intraband contribu- 
tion to <t(w) is reduced as x is reduced. The sharp structures 
at small u> are due to phonons. Interband transitions become 
important for huj > 1 eV. 

To use the f-sum rule, we first make a simple estimate 
of the hopping energy. Due to the doping with x holes 
per site, a given site is occupied by a hole or an electron 
with the probabilities x and (1 — x), respectively. In the 
t — J model, a hole can hop to a neighboring site only if 
the site is occupied by an electron. This suggests a hop- 
ping probability of the order of x(l — x). Multiplying by 
the number of nearest neighbors (four) and the hopping 
integral t gives a simple estimate of the hopping energy 
(Calandra and Gunnarsson, 2003) 



E K = -4tx{l - x)N, 



(16) 



where ./V is the number of sites. An exact diagonalization 
calculation for the t— J model, gives a similar dependence 
on x but a somewhat smaller prefactor of about 3.4 (Ca- 
landra and Gunnarsson, 2003). 

Typically, x -C 1. From Eq. (16) it then follows that 
the right hand side of the f-sum rule (Eq. (6)) should be 



approximately proportional to x, as is also found exper- 
imentally (Yamada et al, 1998). The constant of pro- 
portionality, extracted from Uchida et al. (1991), agrees 
with the calculation to within about 20-30 % (Calandra 
and Gunnarsson, 2003). 

Fig. 12 shows <t(uj) of La2_ a; Sr :! ;Cu04 for different T 
and x (Takenaka et al., 2002). The figure illustrates how 
the peak at u> = gradually disappears as T is increased 
and that for large T, apart from phonon structures, (t(lo) 
behave roughly as assumed in Fig. 10. A comparison of 
the curves for x = 0.06 and x = 0.10, shows that cr(w) 
is much smaller for x — 0.06 in the energy range ui < 1 
eV, where the intraband contributions dominate. This 
illustrates the reduction of the intraband contribution to 
ct(cj) when x is reduced, as assumed above. Similar re- 
sults were found for La2- x (Ca,Sr) x CaCu06+,5 by Wang 
et al (2002). 

We now analyze the resistivity in terms of the f-sum 
rule (Eqs. (6, 8)). For conduction in the ab-plane of 
a cuprate, the factor 1/6 in Eqs. (6, 8) should be re- 
placed by a factor 1/4, since the the system is quasi- two- 
dimensional. Insertion of Eq. (16), but with the prefactor 
3.4, in Eq. (8) gives 



0.4 



Psat — 



0.07c 



x(l — x) x(l — x) 



mficm, 



(17) 



where we have used distance c = 6.4 A between the Cu02 
planes appropriate for La2- x Sr x Cu04. For small x, the 
saturation resistivity in Eq. (17) is very much larger than 
the Ioffe-Regel resistivity, 0.7 mficm. It is also much 
larger than the result in Sec. VI for weakly correlated 
transition metal compounds (Eq. (11)) 



Psat 



O.Ud 
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N~d 



mficm, 



(18) 



where d ~ 3 is expressed in A. This is partly due to 
the degeneracy being just Nd = 1 for the t — J model 
but Nd = 5 for the TM model. Furthermore, the strong 
correlation drastically reduces the hopping energy in the 
cuprates, which gives the factor x(l — x) in Eq. (17). 
Finally, the large separation of the Cu02 planes also in- 
creases p. 

Fig. 13 compares the experimental results of Takagi 
(1992) for L^-rcSr^CuO^ with the saturation resistiv- 
ity in Eq. (17). Due to the factor 1/x in Eq. (17), all 
resistivities are multiplied by the doping x. For small 
x, p sat is then a constant. The experimental resistiv- 
ity is smaller than the predicted saturation resistivity. 
The same conclusion is obtained for other high-T c su- 
perconductors (Calandra and Gunnarsson, 2003). The 
experimental data do therefore not demonstrate absence 
of saturation. On the contrary, for La2- a; Sr a ;Cu04 with 
x = 0.04 and x = 0.07 there are signs of saturation 
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where saturation is expected to occur. It is interest- 
ing that the curves xp{T) for x = 0.04 and x = 0.07 
fall almost on top of each other, giving a dependence 
p ~ 1/x, as suggested by the arguments above. For small 
dopings, saturation has also been reported (Wang et al, 
1996ab) for I^S^Cai-^YzC^Os+j, at similar values as 
for La 2 - a; Sr ;c Cu04 (see Fig. 2). 
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FIG. 13. Resistivity multiplied by the doping x for 
La 2 -*Sr x Cu0 4 (Takagi et al, 1992) for x = 0.04 (full curve), 
x = 0.07 (broken curve), x = 0.15 (dotted curve) and x = 0.34 
(chain curve). The horizontal line shows the saturation resis- 
tivity (Eq. (17)), xp sa , t , in the limit of small x. The figure 
illustrates that the saturation resistivity is not exceeded by 
the experimental data and that there are signs of saturation 
for small x roughly where saturation is expected. 



Above we gave an upper limit to the resistivity with- 
out specifying a scattering mechanism. For x = 0.15, the 
resistivity in the t — J model is of the right order of mag- 
nitude, suggesting that electron-electron scattering can 
explain the resistivity, while for x = 0.07 and x = 0.04 
the resistivity is much too small compared with exper- 
iment (Calandra and Gunnarsson, 2003). This suggests 
that for x = 0.04 and x = 0.07, there is an additional 
important scattering mechanism, beyond the electron- 
electron scattering mechanism in the t — J model, or 
that our small cluster (16 sites) cannot describe some 
important scattering mechanism. This mechanism ap- 
parently increases the resistivity so rapidly for small T 
that the upper limit (17) is approached and the resistiv- 
ity shows sign of saturation. The estimate of the upper 
limit (17) is nevertheless correct, unless the additional 
scattering mechanism appreciably influences the hopping 
energy. The t — J model alone only gives a rather weak 
saturation and only at larger T than seen in Fig. 13. 



VIII. VIOLATION OF THE IOFFE-REGEL CONDITION IN 
OTHER COMPOUNDS 



Figs. 3 and 6 show that the Ioffe-Regel condition is vio- 
lated for several other compounds. We first discuss CrC>2, 
which is a ferromagnet below Tc — 390 K. According to 
band structure calculations, three bands of mainly Cr t 2g 
spin up character are occupied by two electrons (Lewis 
et al, 1997; Mazin et al., 1999). We therefore consider a 
three-fold degenerate band with the width W = 2.3 eV 
(Lewis et al, 1997) and without spin degeneracy. Fig. 
9 compares the experimental cr(u) at T = 300 K (Sin- 
gley et al., 1999) with the large T theoretical model in 
Fig. 10. The theoretical area, ti 2 J a(u)dw = 2.7 (eV) 2 
(for units see Appendix A) has been obtained from the f- 
sum rule (Eq. (6)), assuming a semi-elliptical band with 
W = 2.3 cV. We estimate the experimental t2 S contri- 
bution to the f-sum rule to be 2.2 (eV) 2 , i.e., somewhat 
smaller. As a result, our estimate of the saturation resis- 
tivity (p sat = 0.4 mficm) is somewhat smaller than the 
experimental result (~ 0.6 — 0.7 mSlcm). The reason is 
probably both that the area of cr(u>) is somewhat overes- 
timated and that the theoretical shape of <j{u>) does not 
take into account a strong dip in N{e) at Ep for CrC>2. 
These considerations suggest that theory gives a quali- 
tatively correct description of saturation for CrC>2. The 
strong increase in a(w) for cj <~ 1.5 cV is due to interband 
transitions, which are neglected in the model. 

This treatment assumed Cr02 to be ferromagnetic. If 
we had taken into account that Cr02 is paramagnetic for 
T > Tc, the theoretical saturation resistivity would have 
dropped by a factor of two, due to an increase of \Ek\- 
Experimentally, however, the resistivity changes little at 
T = Tc, suggesting that correlation effects reduce \Ek\ 
above Tc in a way similar to that of magnetic effects 
below T c . Actually, the Hubbard U of Cr0 2 has been 
estimated to U w 3 eV (Korotin et al., 1998), comparable 
to the t2 g band width, and suggesting appreciable corre- 
lation effects. There are rather strong magnetic fluctu- 
ations in several of the systems violating the Ioffe-Regel 
condition. The arguments above suggest that this may 
reduce \Ek\ and increase the saturation resistivity. 

Using similar assumptions (correlation reduces \Ek \ by 
a factor of two in the paramagnetic state) as for Cr0 2 , 
we predict a similar saturation resistivity for VO2 and 
a somewhat larger value (0.6 nrslcm) for CaRu0 3 and 
SrRuC>3. As for Cr02, more experimental data is needed 
to decide whether this agrees qualitatively with exper- 
iment. It is not clear if this framework (strong cor- 
relation effects reducing \Ek\) can be used to under- 
stand Sr2RuC>4, which shows no signs of saturation up 
to T = 1300 K, and La 4 Ru 6 0i9, which shows signs of 
saturation, but at a very large resistivity. 
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IX. VERY LARGE T BEHAVIOR. C 60 COMPOUNDS 

To derive the Ioffe-Regel condition, we assumed that 
ii) T <C W, leading to a T independent upper limit p sat to 
the resistivity. We now consider T to be so large that the 
phonons substantially change W. Such effects are usually 
not very important, but the C@o model is an important 
exception, due to its small band width. The T depen- 
dence of p sat can then be so strong that the concept of 
resistivity saturation becomes meaningless. In the follow- 
ing we neglect the electron-electron interaction, although 
it could play a substantial role for alkali-doped fullcridcs. 
We treat the phonons in the semiclassical approximation 
(Sec. IV). 

For very large T, there is a rather trivial T depen- 
dence due to the Fermi temperature, Tp, entering in the 
Fermi- functions of Eq. (5). Let us consider the resistivity 
due to static disorder. Expanding the Fermi functions in 
l/T, we obtain <r(0) ~ l/T and p{T) ~ T, although the 
scattering mechanism itself is T- independent. A similar 
dependence also enters for the electron-phonon scatter- 
ing, which tends to mask interesting differences between 
couplings to the level energies (LE) and hopping integrals 
(HI). We therefore study the case when the Fermi func- 
tions are replaced by 0-functions in Eq. (5) (Tp = 0). 
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FIG. 14. Resistivity for a model of ordered Ceo molecules, 
considering couplings to the level energies (full line, LE 
coupl.) and to the hopping integrals (broken line, HI cou- 
pling) in the semiclassical treatment of the phonons. Results 
are also shown for Tf = 0, replacing the Fermi functions in 
Eq. (5) by O-functions. The figure shows that there is a large 
difference between LE and HI coupling (after Calandra and 
Gunnarsson, 2002). 

We consider the Cgo model (see Sec. Ill) with either 
LE or HI coupling (Calandra and Gunnarsson, 2002), us- 
ing the same coupling A in both cases (Fig. 14). While 
the resistivity shows no signs of saturation for the LE 
coupling (full curve), the model with HI coupling shows 
a weak saturation (broken curve), i.e., a moderate re- 
duction of the slope for T > 0.03 eV. This saturation be- 



comes much more pronounced for Tp = (dotted curve). 
For the TM model, the slope of p(T) is reduced at high T 
for both HI and LE coupling, but the reduction is more 
pronounced for HI coupling. 

To apply the analysis of Sec. V, we calculate W(T) 
and E K (T). W(T) increases with T for both LE and 
HI coupling. In the LE case this happens because the 
fluctuations in the level positions increase and in the HI 
case because the average of the square of the hopping 
integrals increases. Since the average phonon amplitude 
squared is (x 2 ) ~ T, the second moment S*2(T) of the 
density of states goes as ^(T) = 52 (0) + aT, where a is 
some constant. This leads to (Calandra and Gunnarsson, 
2002) 

w ^ = w ^f+^W§^y (19) 

The T dependence comes in the form \T/W, where 
the large prefactor c = 12 is due to the large ratio 
W 2 1 'Si ~ 12. Using the considerations of Sec. VI, it fol- 
lows that for the HI coupling Ek varies in a similar way, 
and the T dependences of Ek and W essentially cancel 
(Eq. (8)). This is confirmed by Fig. 14 for the C 60 model 
(lowest HI curve) , while the cancellation is less complete 
for the TM model due to a stronger T dependence of 
a (Eq. (10)) and 7 (Eq. (8)). With LE coupling, on 
the other hand, Ek is reduced as T is increased. As the 
phonons are excited, the level energies on different sites 
become different and hopping becomes more difficult. 
For the C 60 model \E K \ ~ 1/^/1 + c\k B T/W(T = 0) 
(Calandra and Gunnarsson, 2002). For the LE coupling, 
the hopping energy and the band width work together, 
and the upper limit to the resistivity takes the form 

^-^ 1+cX wwh) ) m0cm - (20) 

In A3C60 (A= K, Rb) there is orientational disorder, 
i.e., the Cgo molecules more or less randomly take one 
out of two preferred orientations (Stephens et al, 1991). 
As a result, the Drude peak is essentially gone already 
for T = 0. The resistivity can therefore be considered 
as "saturated" already at T = 0. The upper limit for 
the resistivity (Eq. (20)) has, however, such a strong T 
dependence that the term "saturation" becomes mean- 
ingless in this case (see Fig. 7). The reason for this 
strong T dependence is both the small band width and 
the LE coupling to the intramolecular phonons. The T 
dependence of the TM model is weaker, due to the larger 
band width and to the HI coupling. Even in this case, 
however, p(T) does not become a constant because of the 
T dependence of a and 7. 



12 



X. ANDERSON METAL-INSULATOR TRANSITION AND 
MOTT'S MINIMUM CONDUCTIVITY. 

In the semiclassical treatment of the phonons, the 
phonons cause a static disorder, and the problem is there- 
fore related to conduction in disordered system. Thus 
the LE and HI couplings correspond to diagonal and 
off-diagonal disorder, respectively. While the disordered 
systems are usually studied for small T, we are here in- 
terested in the large T behavior. In the semiclassical 
treatment of the phonons, however, apart from causing 
disorder, T only enters via the Fermi-functions (Eq. (5)), 
which is not important for the qualitative behavior. 

Diagonal disorder can lead to an Anderson metal- 
insulator transition at T = (Lee and Ramakrishnan, 
1985). For off-diagonal disorder, however, Antoniou and 
Economou (1977) found that there is no metal-insulator 
transition if the Fermi energy is located in some finite 
region around the middle of the band. The semiclassical 
calculations agree with these results, i.e., localization is 
found for LE but not for HI coupling as T is increased. 

In the QMC calculation there is no sign of localiza- 
tion for LE coupling, just a lack of saturation. This is 
natural. Localization depends sensitively on the phase 
factors, which are not destroyed by the clastic scatter- 
ing in a disordered system. The phase information is, 
however, lost in the inelastic scattering by phonons at 
finite T, and localization is not expected (Lee and Ra- 
makrishnan, 1985). These effects is properly included in 
the QMC but not in the semiclassical treatment of the 
phonons, and therfore localization shows up in the semi- 
classical (Calandra and Gunnarsson, 2002) but not in the 
QMC treatment (Gunnarsson and Han, 2000). 

Mott (1974) has argued that as the disorder increases, 
there is a discontinuous transition from a metal to an in- 
sulator at T = 0. He therefore introduced the concept of 
the minimum conductivity 

a min = 0.026^-, (21) 

where d is the nearest neighbor atomic distance. Later 
work showed that the transition from a metal to an insu- 
lator actually is continuous, but that <r m i n is still relevant 
for low but nonzero temperatures (Lee and Ramakrish- 
nan, 1985). We therefore make a comparison of cr m i n to 
the resistivity in the TM and Ceo models. Converting 
Eq. (21) to a resistivity, we obtain 

Pmax = L6rf mOcm, (22) 

where d is measured in A. Based on experiment, Mott 
deduced a somewhat larger minimum conductivity for 
systems containing transition metal atoms, resulting in 
the maximum resistivity 

Pmax = 1 mficm. (23) 



Mott derived his result for diagonal disorder. His re- 
sult can most naturally be compared with our saturation 
resistivity for HI coupling (off-diagonal disorder), since 
saturation is most pronounced in this case. The resis- 
tivity /9 max is much larger than the saturation resistivity 
obtained above (Eq. (11, 18)) for the TM model with a 
five- fold degenerate orbital (Nd = 5). For a fee lattice 
and a half-filled semi-elliptical band it takes the form 

Asat = -jvT" mStem, (24) 

which is of the order of 0.1-0.2 milcm. The corresponding 
conductivity is substantially larger than Mott's minimum 
conductivity. 

XI. CONCLUSIONS 

We have reviewed experiments showing resistivity sat- 
uration, i.e., p(T) growing more slowly than p(T) <~ T for 
large T. Resistivity saturation is found for several classes 
of metals with large resistivities, in particular for many 
transition metal compounds. Saturation often happens 
in such a way that the Ioffe-Regel condition, I > d, re- 
mains fulfilled. Over the last 15 years, however, a number 
of metals have been found for which the Ioffe-Regel con- 
dition is violated. Some of these metals probably show 
resistivity saturation, but at a much larger value than 
the Ioffe-Regel resistivity, while in other cases little or no 
sign of saturation is seen. 

We have reviewed early theories, presented at a time 
when saturation appeared to be universal. Several of the 
theories derived saturation. These theories emphasized 
different mechanisms for saturation, and no consensus 
was reached about which mechanism dominates. 

Here we argue that it is useful to study the problem 
using the f-sum rule. By assuming that the (Drude) peak 
at u — is gone and that only incoherent contributions 
are left, we obtain an approximate upper limit, p sa t, to 
the resistivity, which usually has a weak T dependence. 
Saturation then happens if the resistivity initially grows 
so rapidly that p sa t is reached for small values of T and 
if p sat has a weak T dependence. 

We have considered three models with qualitatively 
different behavior: 1) A model of weakly correlated 
transition metal compounds, which shows saturation in 
agreement with the Ioffe-Regel condition, 2) a model of 
strongly correlated high-T c cuprates, which can give sat- 
uration but at much larger values than the Ioffe-Regel re- 
sistivity, and 3) a model of alkali-doped Ceo compounds, 
which shows no saturation. 

To derive the Ioffe-Regel condition for model 1) we 
assumed i) noninteraction electrons and ii) T -C W. As- 
sumption i) is violated for model 2) and assumption ii) 
for model 3). The type of electron-phonon coupling, i.e., 
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coupling to the level energies or the hopping integrals, is 
also important when comparing models 1) and 3). 

We have focused on work where cither the electron- 
phonon or the electron-electron interaction was treated 
alone. In many cases this may be an oversimplification. 
For instance, the electron-electron interaction is believed 
to be important for alkali-doped fullcridcs (Gunnarsson, 
1997), although only the clcctron-phonon interaction was 
considered here. Generally, the electron-electron interac- 
tion should reduce the magnitude of the hopping energy 
and tend to increase the resistivity, as was found for the 
high-T c cuprates. The interplay between the electron- 
electron and electron-phonon interactions may, however, 
be more intricate, as found for metal-insulator transi- 
tions (Han et al, 2000) and superconductivity (Han et 
al, 2003) in alkali-doped fullcrides. 



APPENDIX A: MEAN-FREE PATH 

Since the mean-free path and the Ioffe-Regel resistivity 
play an important role in the discussion, we give details 
of how these quantities were obtained. The conductivity 
tensor can be written as (Ashcroft and Mermin, 1976) 



^ 2 E 



d 3 k 



47r3 . T! ,(kK(kK(k)[-g] £=: , (k , ( U) 



where / is the fermi function, r„(k) is the relaxation 
time and v„(k) is the velocity for a state with band in- 
dex v and wave vector k. We assume a three-dimensional 
isotropic system and a spherical Fermi surface with one 
sheet. We furthermore assume that r is independent of 
v and k. Together with vp = hkp/m, where vp is the 
Fermi velocity, and I = tvf, this leads to Eq. (1). As- 
suming that there are M sheets, we find a moderate re- 
duction of the apparent mean free path from Eq. (1) 
by a factor M 1 / 3 . For a quasi- two-dimensional system, 
like the high-T c cuprates, we instead assume a cylindrical 
Fermi surface with the height 2ir/c and radius kp, where 
c is the average separation of the CuC-2 planes. Based 
on photoemission results (Ino et al, 1999; Yoshida et 
al., 2001), we assume a "large" Fermi surface containing 
roughly one electron or hole (more precisely l±x carri- 
ers, where x is the doping). This leads to kp = \/2~7T/a, 
where a is the lattice parameter of the Cu02 plane. Then 



2irhc 



P2d = -o 



e 2 kpl 



(A2) 



Assuming the Ioffe-Regel condition I — a, we find 

p^Jf* = 0.055(c/a ) mficm, (A3) 

where ao = 0.529 A is the Bohr radius and wc have 
used the conversion Hao/e 2 = 0.022 mOcm. Assuming 
c = 6.4 A, appropriate for La2- x Sr x Cu04, we obtain 



P2d = 0-^ m ^ cm - If we instead had assumed a "small" 
Fermi surface, containing x carriers, the resistivity would 
have been 0.7 '/y/x. Even in this case, the experimental 
resistivity of L^-zSr^CuO^ in Fig. 2 exceeds the Ioffe- 
Rcgcl resistivity for small x. 

For A3C60 (A= K, Rb), the Ioffe-Regel resistivity was 
calculated by assuming that I is the separation of two Ceo 
molecules. Scattering inside the molecule is not possible 
at small and intermediate T, since this would involve 
scattering into states that are at least 10000 K higher 
in energy. Therefore the intermolecular separation is the 
appropriate length scale. 
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